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Abstract 

We studied the temperature and the angular dependences of the upper critical field (H C 2(T,0)) 
of Mgi_ x Al x B2 single crystals (x = 0.12 and 0.21) and compared with the dirty-limit two-gap 
theory. We found that H C 2(T,6ys were well described in a unified way by this theory. The 
obtained values of the parameters indicated that as the Al concentration was increased, anisotropic 
impurity scattering increased, making the a bands less anisotropic. Accordingly, the temperature 
dependence of the anisotropy ratio of (7//) systematically decreased, and for x = 0.21, 7^ was 
nearly constant. Our results imply that Mgi_a;Al x .B2 single crystals are in dirty-limit and that 
two-gap nature survives until x = 0.21. 

PACS numbers: 74.70.Ad, 74.62.Dh, 74.62.Bf 
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It is now well established that MgB 2 is a two-gap superconductor with two distinct energy 
gaps: a large gap originatingfrom two-dimensional a bands and a small gap originating from 
three-dimensional ir bands 0, [J. One of the main consequences of the two-gap nature 
is the strong temperature dependence of the H c2 (T) anisotropy, 7^ = H^J H° 2 |4j, which 
is not expected based on the single-gap Ginzburg-Landau theory. Theoretical calculations 
show that the strong temperature dependence of •jh arises from the fact that the anisotropic 
a bands dominate jh at low temperatures while the 7r bands gradually become important 
at temperatures near T c . The above anomalous behavior of 7(T) for MgB 2 single 



crystals was confirmed by using magnetization measurements 

When impurity scattering is increased, the above-mentioned behaviors of H c2 are modi- 
fied. Gurevich Jf]], and Golubov and Koshelev []| formulated the dirty-limit two-gap theory 
for H C 2 by using the quasiclassical Usadel equations. According to this theory, the shape 
of the H c2 {T) curve essentially depends on the diffusivities of the a and the ir bands. For 
T ph T c , H C 2(T) is determined by a maximum diffusivity (cleaner bands) between D a and D n 
while H C 2(0) is controlled by a minimum diffusivity (dirtier bands). When the a bands are 
dirtier, an upward curvature should appear near T c , and 7^ should decrease with temper- 
ature. In contrast, when the it bands are dirtier, a huge increase in H c2 (T) should appear 
at low temperatures without an upward curvature near T c , and 7# should increase with 
temperature. 

Impurity scattering also changes H C 2(8). H C 2(8) was predicted to deviate from the an- 
gular dependence of the anisotropic one-gap Ginzburg-Landau (GL) theory, especially near 
the middle-angle region. This deviation should be most pronounced at T/T c pa 0.95 when 
the parameters supplied by band-structure calculations are used (?[ 13 • Even though these 
predictions were quantitatively compared with H C 2(8) for MgB 2 single crystals and reason- 
able consistency was observed [10], the problem of whether the dirty-limit theory could be 
applied to clean MgB 2 single crystals still remained. In this sense, the dirty-limit theory has 
not yet been verified unambiguously for single crystals in the dirty limit, especially for the 
orientational dependence of H c2 . 

In this paper, we report the effect of Al doping, as deduced from resistance measurements 
at various angles 9 between H and the c-axis, on if c2 (T, 9) of Al-doped MgB 2 single crystals. 
This directional study of the resistance was possible due to success in growing flat and 
regular-shaped Mgi_ x Al :r B 2 single crystals with values of x up to 0.21 and with T c = 25.5 K. 
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We found that two-gap superconductivity in MgB2 was drastically affected by the Al doping 
and that key features predicted by the dirty-limit two-gap theory were observed. Our main 
observations are the following: (1) As the Al concentration increases, the residual resistivity 
(po) greatly increases, implying that Al substitution enhances impurity scattering and that 
the Al-doped samples are in the dirty-region. (2) H c2 {T) can be consistently explained 
within the dirty-limit two-gap theory up to x = 0.21, even though if C 2(0) decreases with 
Al concentration. (3) The 7#(T) systematically decreases and for x = 0.21, 7# is virtually 
temperature-independent. (4) The H c2 {9) for x = 0.12 showed a clear deviation from the 
behavior predicted by the anisotropic GL theory, which is a strong indication of the two-gap 
nature in MgB 2 . However, for x = 0.21, this deviation became very small. The values of the 
obtained parameters suggest that impurity scattering is enhanced in the ir bands, especially 
along the c direction and that the anisotropy of the a bands is significantly reduced. 

Mg] tA LBo single crystals with x = 0.12 and 0.21 were grown under high-pressure con- 
ditions | 111 . Il2j ] and were characterized and patterned as in 

HQ. 

Two sets of samples 

with clean, shiny surface were investigated for each Al concentration. For the resistance 
measurements, well-shaped single crystals with both sides flat were selected from numerous 
samples. The temperature and the angular dependences of the resistance were measured 
from to 9 T by using the AC transport option in a PPMS Quantum Design system. 

Figure [U shows the resistivity p of the Mg 1 _ a ,Al a ,B 2 single crystals (x = 0, 0.12, and 0.21) 
as a function of temperature. As the Al concentration increases, T c decreases. The T c 's 
are 30.8 K and 25.5 K for x = 0.12 and x = 0.21, respectively. The data for x = were 
taken from Ref. |4| and T c of this sample was around 37 K. Previously, for MgB2 single 
crystals, the resistance was reported to follow the Bloch-Griineisen (BG) formula with a 
Debye temperature of 0£> ~ 1100 K |ll(. This implied that the normal-state transport 
properties were well described by an electron-phonon interaction without considering an 
electron-electron interaction. To check whether this is the case in Al-doped single crystals, 
we fitted the p{T) data with the BG formula, where fitting paramters are 0d and residual 
resistivity po- The solid lines in the figure are the BG theoretical curves and describes the 
p{T) data well. The value of G>£> in Al-doped single crystals is ~ 1000 K, which is similar 
to that of MgB 2 single crystals, po increases monotonically with doping, and the fitted 
values of po are 1-63, 21.4, and 32.2 pfl cm for x = 0.0, 0.12, and 0.21, respectively. The 
inset of Fig. ^ shows the normalized low-field magnetization for zero-field-cooled state of 
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Mgi_ :r Al x B2 single crystals extracted from the same batch of single crystals as was used 
for the resistivity measurements. The T c 's determined from the resistivity and from the 
low-field magnetization were virtually the same. 

Figure 121 (a) and (b) show, as an example, the temperature dependences of the resistances 
of the x = 0.12 sample for H \\ c and || ab, respectively. As with MgB 2 single crystals, this 
sample shows surface superconductivity: As the temperature decreases, the resistance first 
decreases linearly and then suddenly drops to zero. In the region of linear decrease, the 
resistance depends on the applied current, and a higher current induces a higher resistance. 
The drop in the resistance indicates the onset of bulk superconductivity. Particularly at high 
currents (1 = 3 mA) for H \\ c, a peak, which is absent at low currents (1 = 1 mA), appears. 
The current dependence of this peak suggests that it is due to the peak effect, observed in 
MgB2 single crystals Q|. The upper critical fields can be determined unambiguously as 
the points where the resistance drops to zero in the curves for I = 1mA. Those points are 
indicated by the arrows. 

In Fig. Efa), H c c2 (T) and H${T) for x = 0.12 and 0.21 are plotted, where H c c2 (T) and 
H^(T) are if C 2(T)'s for H \\ c and for H \\ ab, respectively. For comparison, we also insert 
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H c2 {T) for x = 0.0, which was taken from Ref. |4|. Interestingly, both H^ 2 (T) and H^ 2 (T) 
decrease with increasing Al doping. As a result, the extrapolated H^ 2 (0) and H^ 2 (0) are 
reduced. While the decrease in H® 2 (0) is consistent with the results for polycrystalline 
samples, the decrease in H^ 2 (0) is not. In a study by Angst et al, a small increase in H^ 2 (0) 
was observed at an Al doping of 10 % [14| . By comparing H C 2(0) in both Al- and C-doped 
MgB2, they concluded that in Al-doped samples, the shift in the Fermi level was dominant 
in determining H c2 (T) while in C-doped samples, disorder played a major role. However, 
in light of the huge increase in po, the effects of disorder are not negligible and should be 
taken into account. Another clue to the degree of dirtiness in Mgi_2.ALj.B2 single crystals 
for x = 0.12 and x = 0.21 can be found in the shape of H£ 2 (T) near T c . While MgB2 single 
crystals show a linear decrease in H^ 2 {T) near T c , close inspection reveals that an upward 
curvature gradually appears with Al doping. This becomes even clearer if H^ 2 (T) for x = 
is compared with that for x = 0.21. The upward curvature is consistent with the two-gap 
dirty-limit theory. 

Since the variations in H c2 (T) with Al doping appear to agree well with the two-gju; 
theory, we quantitatively analyzed our H C 2(T) data by using the dirty-limit theory P 
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For x = 0, the dirty limit model may be inappropriate because pure MgB2 crystals are 
considered to be in the clean limit [l5| . If interband impurity scattering is assumed to be 
zero, H c2 (T) for H \\ c is given by 



a [In t + U{h)} [In t + U(r}h)] + a 2 [In t + U(r}h)] 

+ ai [\nt + U(h)] = 0, (1) 

where t = T/T c , U{x) = \&(l/2 + x) — ^(x), ^/(x) is the Euler digamma function, h = 
H c2 D® b /2(f)oT, is the magnetic flux quantum, r\ = D® b /D% b , D*\ is the in-plane electron 
diffusivity of the a and the 7r bands, and ao,i,2 are constants derived from the electron- 
phonon coupling constants (A^J and the Coulomb pseudopotentials (/i m „). The precise 
definitions of 00,1,2 can be found in Ref. 6. For H || ab, the in-plane diffusivities in Eq. Q 
can be replaced by [D^D^] 1 ^ 2 , where T)% m are the out-of-plane electron diffusivities of the 
cr and the 7r bands, respectively. Equation [T] can be generalized to the anisotropic case of 
an inclined field by replacing the diffusivities with the angle-dependent diffusivities D a (6) 
and D T (9) for both bands, where D a ^{6) = [(£>^J 2 cos 2 9 + Df^D^ sin 2 9] 1 / 2 . For the four 
input parameters \ mn = — /i mn at each Al doping level, which reflects the change in the 
electronic structure by electron doping, we used the values determined from first-principle 
calculations [lq . and we obtained the numerical value of the diffusivity for each band. 

In our samples, the interband impurity scattering is believed not to be significant or, if any, 
to be negligible to the first approximation. This is because the interband impurity scattering 
was predicted to eliminate the distinction of each superconducting gap, destructing two- 



gap features Therefore, the upward curvature, which is the hallmark of the two- 

gap superconductivity would have not been observed, if interband impurity scattering were 
significant. 

The solid lines in Fig. Efa) present the theoretical two-gap dirty-limit curves of H c2 (T) 
for x = 0.12 and 0.21. The optimized values of D^ ,c , D® b ' c , and H^iO) from the fits are 
summarized in Table UJ The upward curvature observed near T c for x = 0.12 and 0.21, which 
is typical when a bands are dirtier than ir bands p, may indicate dirtier a bands. If the n 
bands are dirtier than the a bands, the upward curvature near T c should disappear; instead, 
a huge increase in H c2 {T) should appear at low temperatures. The dashed line for x = is 
a guide to eyes. 

Quantitatively, the values of D^ b,c and D® b,c really prove dirty a bands (D^ b,c D^ b,c ), 



which is consistent with the shape of H C 2(T). Dirty a bands were also observed in Al- 
doped MgB2 polycrystalline samples 3|. The electron diffusivity along the c direction in 
the 7r bands is noted to decrease with Al doping while that in the ab plane virtually does 
not change. This originates from pronounced impurity scattering in the n bands as the Al 
concentration is increased. The pronounced impurity scattering, however, is not isotropic as 
is normally assumed. Along the c direction, impurity scattering is more enhanced than in 
the ab plane. Similarly, Al doping influences impurity scattering in the a bands. In this case, 
the electron diffusivity along the c direction increases with Al doping while that in ab plane 
virtually is unchanged. Consequently, the a bands become more isotropic, which is reflected 
in the ratio D a ^ '/D° and this value decreases as Al content increases. The isotropization of 
the a bands is believed to be due to not only the anisotropic impurity scattering but also 
the change in the electronic structure that Al doping induces. 

The same set of electron diffusivities as in Table H] can explain H c2 (8) for x = 0.12 and 
0.21, as shown in Fig. Efb). The solid lines indicate the theoretical curves calculated from 
the dirty-limit two-gap theory. The dotted lines are the theoretical curves of the one-gap 
GL model. The error bars in this data are comparable to or less than the symbol size. The 
two-gap theory describes the data better than the GL model for x = 0.12. For x = 0.12, a 
small difference between the two-gap theory and the GL model is apparent and as predicted, 
is most pronounced at the middle-angle regions. This is a strong indication of the two-gap 
nature of Al-doped MgB2 single crystals. This behavior is very similar to that of MgB2 single 
crystals, where a deviation from GL behavior was observed to be peaked at T »s 0.8T C . For 
x = 0.21, the difference between the two-gap theory and the anisotropic GL model is very 
tiny, as is the case for the temperatures we investigated. Despite the indistinction between 
the anisotropic GL model and the two-gap theory for this doping, the shapes of the H C 2(T) 
curves and the values of the fitted diffusivities guarantee the existence of two distinct gaps. 
If the sample for x = 0.21 followed the one-gap GL model, the upward curvature would not 
be observed. 

Finally, Jh(T) for x = 0, 0.12, and 0.21, extracted from the H C 2(T,9) data, are plotted 
as functions of the reduced temperature T/T c in the inset of Fig. 0(a). The values of 7^ 
are systematically reduced, and for x = 0.21, 7^ is virtually temperature-independent at 
high temperatures, slightly increasing at low temperatures. The at low temperatures 
significantly changes with Al doping and the 7^'s merge to 2 — 2.5 at T = T c for all doping 
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levels. This behavior is thought to result from the isotropization of the a bands. The 
decreasing tendency of jh with increasing temperature for x = 0.12 and 0.21 is in good 
agreement with the case of dirty a bands, predicted by using the dirty-limit two-gap theory. 

If the effects of impurity scattering can be ignored in x = 0.12 and 0.21 single crystals, 
H c2 will evolve according to changes in the electronic structure and in the lattice constant. 
Among these, the main effect is due to changes in the electronic structure caused by doping 
with electrons, resulting in a shift of Fermi level Ep to higher energies. At moderate doping 
levels, where a rigid band model is valid, an increase in Ep modifies the band-averaged Fermi 
velocities, primarily in the a bands and the 7#(0), which is ^ VF = v a p a jvp a in the clean limit. 
Here, vpy is the in- plane (out-of-plane) Fermi velocity of the cr bands. According to the 
calculation by Putti et al. |l9|, v Fa remains approximately constant while v F b a substantially 
decreases with Al doping for x < 0.3. At doping levels of x — 0.0, 0.12, and 0.21, that 
calculation produced r y VF = 5.6, 5, and 4.2, respectively. The value at x — 0.0 is nearly 
consistent with 7#(0) estimated from the experimental data, as shown in the inset of Fig. 
IHJa). In contrast, the values at x — 0.12 and 0.21 are significantly larger than the estimated 
7i?(0). In fact, the 7#(0)'s at x — 0.12 and 0.21 are better represented by the parameter 
7 CT = a/ D^. b /D^, which contains information on not only the Fermi velocity but also impurity 
scattering. Therefore, as we said before, Mgi-^Al^E^ single crystals (x = 0.12 and 0.21) 



are in the dirty limit with anisotropic impurity scattering. 



his is in sharp contrast to 



□ 



14L Il9j . Those polycrystalline 



the conclusions for Al-doped MgB2 polycrystalline samples 
samples might have less impurities than single crystals, which is very improbable in normal 
situations. It is noted that while 7 CT decreases with Al doping, 7^ increases from 1.1 to 1.8 
In fact, the electron diffusivities are related to the value of resistivity by the relation of 
1/p oc N a D a + A^DttpI, where N a and N n are partial densities of state in o and tc bands, 
respectively. In the present case, since the electron diffusivities in the 7r bands are larger 
than those in the a bands, the electron diffusivities in the 7r bands determine the resistivities 
of our samples and resistivity should increase with x. This tendency holds in our samples. 
We calculated the values of resistivities by using the obtained diffusivity values and the 
partial densities of state calculated by Ummarino et al. [16] and obtained 10 and 12 fiQcm 
for x = 0.12 and x = 0.21, respectively. The discrepancy of the absolute values, especially 
for x = 0.21 might originate from a large error in calculating the resistivity of small-sized 
samples. 
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To summarize, we investigated the effect of Al substitution on H c2 (T,9) of MgB 2 single 
crystals. From an analysis of H c2 (T, 6) within the dirty-limit two-gap theory, we found that 
Al substitution influenced the electronic structure complexly; in the n bands, it increased 
impurity scattering along the c direction while it made the a bands less anisotropic. Accord- 
ingly, 7_h-(T) was systematically decreased and for x = 0.21, 7^ was virtually temperature- 
independent. The isotropization, especially of the a bands, originates not only from in- 
creased anisotropic impurity scattering but also from electron doping. In H c2 (9), we also 
observed a strong indication of the dirty-limit two-gap nature of Al-doped MgB 2 . 
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TABLE I: Al content x, upper critical fields i?" 2 (0), electron diffusivities along the ab plane 
(the c axis) in the a and the it bands, D" b ^ and D^ c \ obtained by fitting the H C 2(T) data to 
the dirty-limit model. 
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FIG. 1: Temperature dependence of the resistivity for Mgi_ x Al 2: B2 single crystals (x = 0.0, 0.12, 
and 0.21). The solid lines are theoretical curves of the BG formula. The inset shows the normalized 
low-field magnetization in the zero-field-cooled state. 



FIG. 3: (a) Temperature dependence of H C 2 for Mgi_ :c Al r B2 single crystals (x = 0.0, 0.12, and 
0.21). Open symbols represent H C 2(T) for H\ c and closed symbols represent H C 2(T) for H \\ ab. 
The data for x = 0.0 were taken from Ref. [4]. The inset shows temperature dependence of 7#. 
The open triangle is ^ VF = v^ a /v c Fa , and the open circle and squre are 70- = \/ D^/D^'s. (b) 
Angular dependence of H C 2. The solid lines are the theoretical curves for the dirty-limit two-gap 
model, and the dotted lines are those for the Ginzburg-Landau theory. 



FIG. 2: Temperature dependence of the resistance for (a) H \\ c and (b) H \\ ab. 
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Figure 3(a) by Kim etal. 
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Figure 3(b) by Kim etal. 
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